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Abstract 

We present the energy spectra of the fluxes of positrons, anti-protons 
and photons generated by Dark Matter annihilations in our galaxy, as 
univocally predicted by the model of Minimal Dark Matter. Due to 
multi-TeV masses and to the Sommerfeld enhancement of the anni- 
hilation cross section, distinctive signals can be generated above the 
co background, even with a modest astrophysical boost factor, in the 

range of energies soon to be explored by cosmic ray experiments. 

co 
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1 Introduction 

oo 

We consider Minimal Dark Matter [1, 2] (MDM), i.e. we assume that the DM is the neutral 
component of one single weak multiplet, that interacts with SM particles only via (broken) 
gauge SM interactions. The assignment of spin and SU(2)l x U(1) y quantum numbers fully 
identifies each different MDM candidate: the full list, together with a short list of the most 
interesting candidates, has been presented in [1, 2]. 

The main virtues of such model (which is not inspired by more ambitious beyond-the-SM 
constructions like super symmetry or extra dimensions) can therefore be summarized in terms 
of economy and predictiveness. The model has no free parameters as all DM couplings are 
predicted by gauge invariance and the DM mass is determined by matching the relic abundance, 
Qdm^ 2 = 0.110±0.005 [3]. A particularly interesting MDM candidate is the fermionic SU(2)l 5- 
plet with hypercharge Y = 0, that is automatically stable on cosmological time-scale thanks to 
the SM gauge and Lorentz symmetries, without having to impose ad-hoc parities (like i?-parity 
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Figure 1: Velocity dependence of Sommerf eld-enhanced MDM annihilation cross sections, for 
the two candidates that we mainly consider. 



or KK-parity). We will restrict our study to three particularly interesting MDM candidates: 
the fermionic quintuplet already mentioned above; the fermionic 3-plet with hypercharge Y = 
(the MDM candidate that has the same quantum numbers of the supersymmetric 'wino'); the 
scalar triplet with Y — 0. 

We here compute the "indirect DM signals", generated by DM DM annihilations into p, 
e + , 7 in our galaxy [4]. Unlike DM candidates motivated by a successful natural solution 
to the hierarchy problem (which should therefore have mass below or around the Z mass), 
MDM predicts specific multi-TeV values for the DM mass Mom, and annihilation cross sections 
enhanced by electroweak Sommerfeld corrections: these two features imply a distinctive DM 
signal that can be tested by running and future experiments like PAMELA [5] and AMS- 
02 [6], dedicated to extending our knowledge of galactic Cosmic Ray (CR) spectra up to higher 
energies. 



2 Energy spectra at production 

MDM annihilates at tree level into W + W~ , and at loop level into 77, 7Z, ZZ. The relative 
cross-sections are significantly affected by non-perturbative Sommerfeld corrections [7] , and we 
use the results of [2]. As a consequence of Sommerfeld corrections, the DM DM annihilation 
cross sections exhibit a quite steep dependence on Mom and can vary by about one order of 
magnitude within the range allowed at 3cr by the cosmological DM abundance as computed 
assuming thermal freeze-out (see e.g. fig.s 2-5 in [2]). For the same reason, the cross section 
o~v also depends on the DM velocity v, reaching a maximal value for v — > 0, as shown in fig. 1. 
The average DM velocity in our galaxy, v ~ 10~ 3 , is however low enough that av is close to its 
maximal value, which we assume. 

For definiteness, in the following we adopt the following best-fit values of Mdm and (av) for 
the MDM candidates that we consider: 
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Table 1: Coefficients for the analytic approximation in eq. (5) to the fragmentation functions. 
for the fermion quintuplet with Y — 0, and 
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for the fermion triplet with Y — 0, and 
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for the scalar triplet with Y = 0. We will not plot predictions to the scalar triplet, that 
can be easily read out from the corresponding predictions for the fermion triplet, taking into 
account that they have a similar mass, and multiplying all rates by a factor of about 16 [2], 
due to a large Sommerfeld enhancement. The other automatically stable MDM candidate, the 
scalar eptaplet, is expected to have a mass Mdm ~ 25TeV, but we cannot reliably predict its 
annihilation cross sections. 

Annihilation cross sections into and ZZ are given by 

a^z = 2<7 77 / tan 2 #w = 6.5cx 77 , o~zz = 0>y/ tan 4 ^ = 10.8cr 77 . (2) 

for all MDM candidates with Y = 0. 

We next need to compute the energy spectra of e + ,p~ , 7 produced by decays of SM vectors. 
Instead of using the results available from the literature, we performed an independent compu- 
tation. Indeed, to our knowledge, previous 'decay' computations do not take into account spin 
correlations of SM vectors V in the intermediate state. We instead compute the full matrix 
element for DM DM — > VV — > 4 fermions. For example, including spin correlations the energy 
spectra of primary positrons directly produced in the annihilation DM DM — > W + W~ — > 4 
fermions, is 
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where x = E e + /Mom and Mdm 3> Mw- The same x-dependence applies to all massless fermions 
and SM vectors, and it arises as follows (see e.g. [8]). SM vectors are produced isotropically in 
the DM DM rest frame, with equal transversally polarized helicities h: both + or both — . This 
is a characteristic of s-wave annihilations, that dominate in the non relativistic limit. A V at 
rest decays into massless fermions with angular distribution 

' /A oc(l + cos#) 2 (4) 
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where cos 9 is the angle between the direction of the fermion and the spin of the vector 
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Figure 2: Energy spectra ofe + ,p,'y produced by non-relativistic DM DM annihilations into SM 
vectors. Only e + have a secondary component (dashed green line shown on the W + W~ plot), 
that dominates at large x ~ 1. 



We generated 4 fermions events with the MadGraph 4.2 [9] MonteCarlo event generator, 
where we extended the SM to incorporate the DM particle and its charged partner. This 
extension incorporates proper tree level interactions among the W and the new particles. The 
interactions of the new fields with the Z boson arises at one loop, and has been taken into 
account effectively adding a tree level vertex between the Z and the neutral DM. Although not 
physical, this vertex is suitable for an easy implementation in the model and at the same time 
provide realistic DM DM — > ZZ collisions. In fact we checked that it leads to the expected 
isotropical production of Z pairs with same helicity, thus proving it to be equivalent to the real 
one loop vertex. Actually, in the case of DM DM — > W + W~ — > 4 fermions, we instead used a 
MonteCarlo routine written by us. 

The decay products of the Ws and the Zs produce QED and QCD radiation. This emission 
is simulated through the parton shower MonteCarlo Pythia 8.1 [10], in which we allowed 
emission from all the final state particles and resonances. This does not take into account the 
emission from intermediate states of the process as simulated at the matrix element level. This 
emission is discussed below. 

The formation of hadronic states out of the shower's products has been simulated with 
Pythia as well. Since we are interested to the observation of stable particles we explicitly 
requested the decay of /i , r and of all unstable mesons and baryons (including n and n). 

The whole showering, hadronizations and decay process results in the production of e + , p~ , 7 
with lower energy x < 1. 

The final energy spectra are plotted in fig. 2 and table 1 provides the numerical coefficients 
c n in the analytic approximation 

(5) 

These results apply in the limit M DM 3> M wz and have no logarithmic dependence on M DM , 
since the virtuality of final state particles is ~ My/,z <C Mdm- 

As Pythia only takes into account brehmstrahlung from the final state fermions, we must 
separately add the photons produced by brehmstrahlung from W and from the fermionic 
charged components DM (relevant at x<l) [11]. Both particles have virtuality ~ Mdm, 
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Figure 3: Photon flux from the galactic center for the isothermal (dashed, J = 13.5) and 
NFW (dotted, J = 1380j DM density profiles. 
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where x = E^/Mbu < 1 — e. 

We next need to consider how 7, e + ,p are produced and propagate in our galaxy. 



3 Photons 

3.1 Astrophysics 

We consider three possible DM halo profiles: cored isothermal [12], the Navarro- Frenk- White 
(NFW) [13] and Moore [14]. In all cases the DM density profile can be parameterized as 



P(r) = P© 
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where r Q = 8.5 kpc is the Earth distance from the galactic center, p & = p(r Q ) is the DM 
density at the Earth position (we assume p Q = 0.3 GeV/cm 3 : values in the range 0.2 — 0.7 are 
considered in the literature [15] and the 7,e + ,p fluxes scale as p 2 ) and the a,f3,j,r s profile 
parameters are: 
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As well known, the NFW and Moore exhibit a cusp at the center of the galaxy. 2 
Photons propagate freely. The differential flux of photons received from a given angular 
direction dQ is 

rf$ 7 1 c pi , dNl t f ds 



dfidE 2 4ttM 2 m ^ d£ ' Anc-of-sight r VP© 

where the adimensional quantity J encodes the astrophysical uncertainty. When observing a 
region with angular size Q the factor J dQ gets replaced by JQ — f fi J dQ. For Q = 1CT 3 
centered around the galactic center one has J = 13.5 for the isothermal DM profile, J = 1380 
for the NFW profile, J = 3830 for the Moore profile. 



3.2 Results 

Fig. 3 shows the predicted energy spectrum of the flux of photons produced by MDM annihila- 
tions. It somewhat differs from the analogous figure in [2] because we more precisely computed 
the photon spectrum. We assumed realistic detector parameters: an energy resolution of 15%, 
and that the region observed has angular size Q = 10~ 3 centered around the galactic center. 
It can be rescaled to any other search strategy. For example, the H.E.S.S. [17] experiment has 
a much better angular resolution, that allows it to resolve the black hole at the center of our 
galaxy. We do not here address which observational strategy maximizes the sensitivity to MDM 
photons: focus on the black hole or subtract it; focus on regions of the galaxy far from the 
center that have less astrophysical 7 sources or on nearby galaxies. These choices only affect JQ 
and not the energy spectrum of MDM photons. The scalar eptaplet with mass Mdm ~ 25 TeV 
provides an energy spectrum that resembles the one emitted by the galactic center. 



4 Positrons 

4.1 Astrophysics 

The positron flux per unit energy from DM annihilations in any point in space and time is 
given by $ e +(i, x, E) = v e +f/4iT (units 1/ GeV • cm 2 • s • sr) where v e + is the positron velocity 
(essentially equal to c in the regimes of our interest) and the positron number density per unit 
energy, f(t,x,E) = dN e +/dE, obeys the diffusion-loss equation: 

f t - K{E) ■ V 2 / - — (b(E)f) = Q (10) 

with diffusion coefficient K{E) = K {E/GeV) s and energy loss coefficient b(E) = E 2 / ( GeV • 
te) with te = 10 16 s. They respectively describe transport through the turbulent magnetic 
fields and energy loss due to synchrotron radiation and inverse Compton scattering on CMB 



2 In various numerical computations, it is convenient to smooth out this behavior adopting the prescription 
discussed in [16]. It simply amounts to replace the divergent profile by a well behaved one below an arbitrarily 
chosen critical radius of r cr i t = 0.5 kpc from the galactic center, while preserving the absolute number of 
annihilations in that region. More precisely, we use 
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Figure 4: Left: The uncertain 'halo function' I(Xd) of eq. (13) that encodes the astrophysics 
of DM DM annihilations into positrons and their propagation up to the Earth. The diffusion 
length is related to energy losses as in eq. (14). Right: The p astrophysical function R(T) of 
eq. (23), computed under different assumptions. In both cases, the dashed (solid) [dotted] bands 
assumes the min (med) [max] propagation configuration of eq. (11) and eq. (22) respectively. 
Each band contains 3 lines, that correspond to the isothermal (red lower lines), NFW (blue 
middle lines) and Moore (green upper lines) DM density profiles. 



photons and on infrared galactic starlight. Eq. (10) is solved in a diffusive region with the shape 
of a solid flat cylinder that sandwiches the galactic plane, with height 2L in the z direction 
and radius R = 20 kpc in the r direction [18]. The location of the solar system corresponds 
to x = (r ,z ) = (8.5 kpc, 0). The boundary conditions impose that the positron density / 
vanishes on the surface of the cylinder, outside of which positrons freely propagate and escape. 
Values of the propagation parameters 5, K and L are deduced from a variety of cosmic ray 
data and modelizations. We adopt the sets discussed in [19]: 

Model 5 Kq in kpc 2 /Myr L in kpc 



min (M2) 0.55 0.00595 1 , , 

med 0.70 0.0112 4 ^ ' 

max (Ml) 0.46 0.0765 15 

Finally, the source term due to DM DM annihilations in each point of the halo with DM 
density p(x) is 

' k 

where k runs over all the channels with positrons in the final state, with the respective thermal 
averaged cross sections av. 

One assumes steady state conditions, so that the first term of eq. (10) vanishes, and the 



7 



Halo model 


Propagation 


a 


a i 


a 2 


a 3 


bi 


b 2 


C\ 


c 2 




/l\ /TO \ 

mm (M2) 


0.500 


0.774 


-0.448 


0.649 


0.096 


192.8 


0.211 


33.88 


NFW 


med 


0.502 


0.621 


0.688 


0.806 


0.891 


0.721 


0.143 


0.071 




max (Ml) 


0.502 


0.756 


1.533 


0.672 


1.205 


0.799 


0.155 


0.067 




mm [nlz ) 


U.oUU 


u. i yi 


n A A Q 
-U.44o 


U.ooo 


U.Uyo 


1 no s 


n 01 1 
U.zll 


66. oo 


Moore 


med 


0.503 


0.826 


0.938 


0.610 


0.912 


0.762 


0.162 


0.055 




max (Ml) 


0.503 


0.889 


1.778 


0.571 


1.230 


0.811 


0.135 


0.061 




min (M2) 


0.500 


0.903 


-0.449 


0.557 


0.096 


192.8 


0.210 


33.91 


isoT 


med 


0.495 


0.629 


0.137 


0.784 


0.766 


0.550 


0.193 


0.296 




max (Ml) 


0.499 


0.695 


0.677 


0.721 


1.092 


0.951 


0.379 


0.231 



Table 2: Fit parameters for the expression in eq. (17) for the halo function I(Xd) that encodes 
the astrophysics of the production density and the propagation of positrons in the galactic halo. 



solution for the positron flux at Earth can be written in a useful semi-analytical form [19, 20]: 
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where B > 1 is an overall boost factor discussed below, Xd{E,E') is the diffusion length from 
energy E' to energy E: 



X 2 D = AKqTe 
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6^1 



5-1 



(14) 



and the adimensional 'halo function' I(Xu) [19] fully encodes the galactic astrophysics and is 
independent on the particle physics model. Its possible shapes are plotted in fig. 4 for the 
set of DM density profiles and positron propagation parameters that we consider. 3 From the 
numerical computation we find that I(Xd) is well reproduced with a naive fit function of the 
form 



H^d) = a o + a-i tanh 
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with £ = log 10 Xjj/kpc an d the coefficients reported in table 2. 

The main features of the halo function can be understood as follows. It is defined such that 
/ ~ 1 at Ad C r Q ,L: all positrons created close enough to the Earth can reach it without 
loosing energy. / can exhibit a peak at ~ r if positrons produced by DM DM annihilations 



3 Formally, one finds that 
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where Jj is the Bcsscl function of the first kind (cylindrical harmonic) of order i, C„ is the n-th zero of the i = 
function and R n ,m corresponds to the Bessel- and Fourier-transform of (p/ Pq) 2 ' 
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Figure 5: Positron fraction, N e +/(N e + + N e -), generated by DM DM annihilations. The red 
(upper) curves refer to the 5-plet MDM candidate (eq.lb). The blue (lower) ones to the 3-plet 
(eq.la). In the left plot we fix the NFW halo profile and vary the e + propagation model. In the 
right plot we fix the med propagation model and vary the DM halo profile. We assumed a boost 
factor B = 10: notice that a signal above the background is present even for B = 1, for the 
5-plet case. The experimental data points are taken from [25, 26, 21, 28]. 



around the galactic center are dominant and reach us after loosing some energy. If instead the 
diffusive region is thin we only receive positrons produced within a region ~ L around the 
Earth: e.g. the dashed lines are for L ~ 1 kpc. 

4.2 Results 

Unlike photons, where we look at the central cuspy region of the galactic DM halo (such that 
the photon flux can be very large but also very uncertain), positrons do not have directionality. 
Especially at energies just below the DM mass M, positrons are dominantly produced close to 
the solar system, so that their flux is less affected by uncertainties in the DM profile. 

However, the DM density in our galaxy might have local clumps that would enhance the 
positron flux by an unknown 'boost factor' B > 1. We take it as energy independent and 
with a value of B = 10. This is a simplifying (but widely used) assumption. Detailed recent 
studies [22, 23, 24, 21] find that a certain energy dependance can be present, subject to the 
precise choices of the astrophysical parameters. Within the uncertainty, these studies also 
converge towards small values of B (except for extreme scenarios), with B = 10 still allowed. 

The results are shown in term of the energy spectrum of the positron flux at Earth from DM 
DM annihilations, computed for several astrophysical models and compared with the expected 
background. The latter, believed to be mainly due to supernova explosions, is obtained by CR 
simulations [29] and can be parameterized as described in [30] by $g+ g = 4.5 E ' 7 / (1 + 650 E 2:3 + 
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1500 E 42 ) for positron and 

$ bk g = $ bkg, P rim +(|) bk g ,sec = o.l6 J B- L1 /(l+llE a9 +3.2 J E; 2 - 15 )+0.70 J E; a 7(l + 110 J E; L5 +580 J E; 4 - 2 ) 

for electrons, with E always in units of GeV. In fig. 5 we actually plot the positron fraction, 
$ e +/($ e + + $ e -)j as the flux $ e - of cosmic ray electrons provides a convenient normalization, 
and the ratio does not depend on solar activity (see the discussion in the case of anti-protons). 
We see that in the region at E ~ Mdm/3 where the signal/background ratio is maximal, the 
predicted signal is quite distinctive and does not significantly depend on unknown astrophysics, 
being manly generated by prompt positrons created close to the Earth. 

The overall rate is however uncertain because the DM DM annihilation cross-sections vary by 
about one order of magnitude within the narrow range of Mdm that reproduces the cosmological 
DM density and because of the possible enhancement coming from the boost factor B. We here 
assumed the sample values of eq.s (1) and a boost factor: B = 10. With this choice, the excess 
starts to appear just around the maximal energy probed by current experiments, and would 
give a clear signal in the PAMELA experiment [5]. 

A signal is still present even for B — 1; however it appears only at higher energies, around 
the peaks of the signal curves in fig. 5 at E m Mdm/3 ~ TeV. This region of energies will 
be hopefully explored by the future AMS-02 experiment [6]. 4 Finally, we recall that the scalar 
triplet MDM predicts a signal about 16 times higher than the fermion triplet MDM. 



4.3 Synchrotron radiation 

Another possible DM signal is the synchrotron radiation from e produced in DM DM anni- 
hilations. For simplicity, since astrophysics is anyhow significantly uncertain, we neglect the 
time and space dependence of / in the diffusion equation (10): this amounts to assume that e 
are trapped enough in the galaxy bulge that they loose there most of the energy, consequently 
maximizing the synchrotron signal. Solving eq. (10) the e energy spectrum is then given by 

1 f°° 

f{E) = W)J E dE ' Q{E '^ E = xM ™, (18) 
and the energy spectrum of synchrotron radiation is then given by 

HP \/3p 3 B f 1 

Eg. oc yA e -^ / Hx f(E = xM DM )F(r/x 2 ) (19) 
HE y 2ir m e J Q 

where F(x) = x K 5 / 3 (^)H^ ~ x 1 ^ 3 e~ x is the synchrotron function. The adimensional factor 
r = 2m 3 ,E 1 /3eBMp M encodes the dependence on E 7 , on the DM mass M DM and on astrophysics 
trough the uncertain magnetic field B. Numerically r pa 10~ 4 for a magnetic field B = fiG, 
Mdm = 1 TeV and E 1 = 10 GHz. The integral in eq. (19) is easily computed numerically; 
we here just notice that one roughly has f(x) oc 1/x 2 due to the E 2 dependence of b(E) and 
consequently HP y /HE 1 oc E 7 1 ^ 2 . 

WMAP observed of an apparent excess of radiowaves with v ~ 20 GHz from the galactic 
center, that might be due to synchrotron radiation from e produced in DM DM annihila- 
tions [32]. The angular dependence of the signal is precisely measured; however it depends on 

4 The experimental limitation on maximal energies arises because the energy is measured from deflection of 
charged particles in the magnetic field of the spectrometer and above a certain threshold positrons are shadowed 
by the abundant spillover protons [31]. The PAMELA experiment should soon release data about e + and p up 
to about 190 GeV and 270 GeV respectively, while AMS-02 might reach the TeV region. 
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astrophysical issues: the DM density profile and the e propagation model. A cusped halo 
model allows to fit the anomaly. The energy dependence of the signal has not been precisely 
measured and the MDM prediction is compatible with the WMAP haze. Furthermore for the 
MDM values of av and of Mdm, and for reasonable values of the magnetic field B, the intensity 
is comparable with the WMAP haze [32]. 



5 Antiprotons 

5.1 Astrophysics 

The propagation of anti-protons through the galaxy is described by a diffusion equation anal- 
ogous to the one for positrons. Again, the number density of anti-protons per unit energy 
f(t, x, T) = dNp/dT vanishes on the surface of the cylinder at z — ±L and r = R. T = E — m p 
is the p kinetic energy, conveniently used instead of the total energy E (a distinction which will 
not be particularly relevant for our purposes as we look at energies much larger than the proton 
mass m p ). Since m p ^> m e we can neglect the energy loss term, and the diffusion equation for 
/is 

df d 

-± - K(T) ■ V 2 / + — (sign(z) / Konv) =Q-2h 5{z) T ann f (20) 



where: 



The pure diffusion term can again be written as K{T) = K (3 (p/ GeV) 5 , where p = 

(T 2 + 2m p T) 1//2 and (3 = v p /c = (l — m 2 /(T + m p ) 2 ) 1 '' 2 are the antiproton momentum 
and velocity. 5 and K are given in eq. (22). 

The V conv term corresponds to a convective wind, assumed to be constant and directed 
outward from the galactic plane, that tends to push away p with energy T < 10 m p . Its 
value is given in eq. (22). 

The source term Q due to DM DM annihilations has a form fully analogous to eq. (12), 
with E now formally replaced by T. 

The last term in eq. (20) describes the annihilations of p on interstellar protons in the 
galactic plane (with a thickness of h = 0.1 kpc <C L) with rate r ann = («H+4 2 ' /3 nHe)c r p^ n Wp, 
where nn ~ 1/cm 3 is the hydrogen density, «He ~ 0.07 tt-h is the Helium density (the factor 
4 2//3 accounting for the different geometrical cross section in an effective way) and is 
given by [33, 20] 

ann _ f 661 (1 + 0.0115 T~ - 774 - 0.984 T 0151 ) mbarn, for T < 15.5 GeV . . 
a PP I 36 T^ 5 mbarn, for T> 15.5 GeV ' ^ ' 



- We neglect the effect of "tertiary anti-protons" . This refers to primary p after they have 
undergone non-annihilating interactions on the matter in the galactic disk, losing part 
of their energy. The effect can be included in terms of an absorption term analogous to 
the last term of eq. (20) but proportional to a different <j non - ann ^ anc [ Q f a re-injection 
term Q tert proportional to the integrated cross section over f(T). The full solution of 
the resulting integro-differential equation can be found in [34] . The effect of tertiaries is 
mainly relevant at low energies T < few GeV. 
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The set of propagation parameters in the case for anti-protons that we adopt has been 
deduced in [35] from a variety of cosmic ray data and modelization (see [36]): 



Model 


5 


Kq in kpc 2 /Myr 


L in kpc 


V conv in km/s 


min 


0.85 


0.0016 


1 


13.5 


med 


0.70 


0.0112 


4 


12 


max 


0.46 


0.0765 


15 


5 



(22) 



Assuming steady state conditions the first term in the diffusion equation vanishes, and the 
equation can be solved analytically [37, 38, 36]. In the "no-tertiaries" approximation that we 
adopt, the solution for the antiproton flux at the position of the Earth &p(T, r ) = Vp/(A7x)f 
acquires a simple factorized form (see e.g. [35]) 



B — 

4:71 



Pe 
Mi 



DM 



. 1 (!N k 



dT 



(23) 



where B is the boost factor. The k index runs over all the annihilation channels with anti- 
protons in the final state, with the respective cross sections; this part contains the particle 
physics input. The function R(T) encodes all the astrophysics and is plotted in fig. 4 for 
various halo and propagation models. 5 From the numerical computation we find that R(T) is 
well reproduced with a fit function of the form 



log 10 \R(T) /Myr] = a + a x r + a 2 r 2 + a 3 r 3 + a 4 r 4 



(26) 



with r = log 10 T I GeV and the coefficients reported in table 3. 

Finally, for completeness we also take into account the average solar modulation effect, 
although it is relevant only for non-relativistic p: the solar wind decreases the kinetic energy 
T and momentum p of charged cosmic rays such that the energy spectrum dQp^/dT^ of anti- 
protons that reach the Earth with energy T and momentum p e is approximatively related to 
their energy spectrum in the interstellar medium, dQp/dT, as [39] 



d$ f 



dTrr 



p% d<5> p 



P 



dT 



T = T (S + \Ze\ 



P 



2m v T + T 



(27) 



The so called Fisk potential 4>f parameterizes in this effective formalism the kinetic energy loss. 
A value of <f>F = 0.5 GV is characteristic of a minimum of the solar cyclic activity, corresponding 
to the period in which most of the observations have been done in the second half of the 90 's 
and are being done now. 



^Formally, it is given by 



oo 

R(T) = J2 Jo (Cn^)exp 



2K(T) 



Vn{L) 
A n smh(S n L/2) 



with 



Vn{Z) 



JKCn)R 2 Jo 



dr r Jo(Cnf I R) / dz exp 



Vconv(Z — z) 

2K(T) 



sinh (S n (Z-z)/2) 



Pq 



(24) 



(25) 



The coefficients A n = 2hT ann + V conv + K(T) S n coth(S n L/2) with S n = (V? onv /K(T) 2 + 4(?JR 2 ) 1/2 encode 
the effects of diffusion. 
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Halo model 


Propagation 




a 


a i 


a 2 


a 3 


a.i 




min 


0.913 


0.601 


-0.309 


-0.036 


0.0122 


NFW 


med 


1.5 


60 


0.517 


-0.293 


-0.0089 


0.0070 




max 


2.740 


-0.127 


-0.113 


0.0169 


-0.0009 




min 


O.J 


94 


0.606 


-0.299 


-0.041 


0.0128 


Moore 


med 


li 


70 


0.553 


-0.289 


-0.0149 


0.0079 




max 


2.5 


10 


-0.119 


-0.117 


0.0181 


-0.0010 




min 


0.927 


0.590 


-0.315 


-0.0319 


0.0115 


isoT 


med 


1.790 


0.399 


-0.315 


0.0162 


0.0031 




max 


2.480 


-0.156 


-0.098 


0.0132 


-0.0005 



Table 3: Fit parameters for the expression in eq. (26) for the propagation function R(T) that 
encodes the astrophysics of the production density and the propagation of antiprotons in the 
galactic halo. 

5.2 Results 

Fig. 6 shows the results for final p flux at earth (at the top of the atmosphere) from DM DM 
annihilations, compared to the background and to the currently available experimental data. 
The background is borrowed from the detailed analysis in [45] , the results of which we find to 
be well reproduced by a fitting function of the form 

log 10 $p kg = -1.64 + 0.07r-r 2 -0.02r 3 + 0.028 r 4 

with r = log 10 T/ GeV. We take for definiteness the flux corresponding to the 'med' propagation 
parameters; see [45] for a complete discussion on the effects of changing that. Particularly 
favorable is the fact that the uncertainty in the estimates of the background is quite narrow 
around 10 — 100 GeV, where results are expected soon. 

The shape of the spectrum appears to be relatively independent from the propagation 
model (fig. 6a) and the halo profile (fig. 6b). Different p propagation models instead change the 
overall signal rate by about one orders of magnitude, consistently with previous results in the 
literature [45, 34, 35]. Different halo profiles with fixed p & make only a difference of a factor of 
a few, which can be interpreted in terms of the fact that the signal is not dominated by the far 
galactic center region, where profiles differ the most. 

As for the case of the positrons, we have plotted the results assuming a modest and energy 
independent boost factor B = 10 (in principle this boost factor and its properties are different 
from those for positrons [23]). In this case the excesses appear in the range of energies soon 
to be explored. Even for a boost factor B = 1 a signal is present above the background for 
most choices of parameters, although it would show at higher energies. Again the scalar triplet 
predicts a signal about 16 times larger than the fermion triplet. 

6 Conclusions 

We computed the indirect detection signatures (fluxes of positrons, anti-protons, 7 and syn- 
chrotron radiation) as predicted by Minimal Dark Matter. We focussed on three particularly 
interesting MDM candidates: the automatically stable ermion 5-plet with hypercharge Y = 0; 
the wino-like fermion 3-plet with Y = 0; its scalar analogous. We fixed the MDM masses and 
annihilation cross sections to the central values predicted in terms of the measured cosmologi- 
cal abundance. These values are listed in eq. (1). Since MDM predicts multi-TeV masses and 
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DM halo model: NFW 



p propagation model: med 




Figure 6: The antiproton flux generated by DM DM annihilations for the case of the fermion 
3-plet (blue, higher) and 5-plet (red, lower) MDM candidates, compared with the astrophysical 
p background (shaded area) and experimental data. In the left plot we fix the NFW halo profile 
and vary the p propagation model. In the right plot we fix the med propagation model and vary 
the DM density profile. The compilation of data point includes results from the BESS [40], 
MASS [41], CAPRICE [42] and AMS-01 [43] experiments, as well as the preliminary results 
from the PAMELA experiment [44]- A boost factor B = 10 is assumed, but a signal is present 
even for B = 1. 



the Sommerfeld electroweak non-perturbative enhancement of the DM DM annihilation cross 
sections into H^+VT - , 77, -yZ, ZZ, the signals for indirect detection turn out to be distinc- 
tive, reaching multi-TeV energies and being above the astrophysical background. The spectral 
shapes are characteristic of DM DM annihilations into SM vectors. 

We recomputed independently most of the ingredients that are necessary for the analysis, 
finding agreement with results in the literature when these are available. The spectra at produc- 
tion were recomputed taking into account spin correlations among SM vectors, by implementing 
the MDM interactions in MadGraph [9] or by a custom-built MC routine and hadronizing 
the resulting MonteCarlo events with Pythia [10]. The propagation diffusion-loss equation for 
positrons and anti-protons was solved following the semi-analytic prescriptions discussed in the 
literature and for a variety of propagation models and halo profiles. We offer plots (fig. 2 and 
4) and simple fit functions for all these ingredients. 

Fig. 3 shows the predicted photon energy spectrum. 

The positron flux is shown in fig. 5. This signal is only very mildly affected by the DM 
density profile, except for multiplicative uncertainties due to the uncertain local DM density 
(we assumed p Q = 0.3GeV/cm 3 ) and due to the boost factor B > 1, taken for simplicity as 
energy independent. On the other hand, e + fluxes somewhat depend on the e + propagation 
model in our galaxy; this uncertainty will be reduced by future measurements of cosmic rays 
and is present only at E <C Mdm- Indeed we do not know if positrons produced around the 
galactic bulk escape from the galaxy or reach us after loosing most of their energy (giving also 
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an interesting synchrotron radiation signal). On the contrary positrons produced in the region 
of the galaxy close to the solar system surely reach us without loosing significant energy: such 
positrons with E < Mdm give a detectable MDM signal even if B = 1. 

The anti-proton flux is shown in fig. 6. Different p propagation models give p fluxes that 
differ by about one order of magnitude, and different DM density profiles give fluxes that differ 
by about a factor 2. A detectable signal is again typical, and its energy spectrum (as predicted 
by MDM) is not significantly affected by astrophysical uncertainties. 

Experimental results are expected soon. 
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Figure 7: Predicions of the Minimal Dark Matter fermion quintuplet compared to preliminary 
PAMELA data. In both figures the boost factor equals B = 3 if (av)ww = 9 • lCT 23 cm 3 /sec 
which corresponds to M = 9.2 TeV, or B = 30 if (o~v)ww — 11 ' 10 -23 cm 3 /sec ; obtained for 
M = 9.6 TeV ; within the range inferred from the measured DM cosmological abundance. 

Addendum: preliminary PAMELA results 

The PAMELA collaboration presented preliminary data about the spectra of anti-proton and 
positron cosmic rays. No anomaly seems present in the anti-proton data, see fig. 7b [1A]. 
Compatibly with measurements from previous experiments, an anomaly possibly due to DM 
annihilations seems present in the data about the e + /(e + + e~) fraction, as indicated by the 
fact that the positron fraction grows at high energy, see fig. 7a [1A]. At energies below 10 
GeV PAMELA measures a flux somewhat lower than previous experiments: this seems due 
to the variation in solar activity [1A]. Previous experimental results have significantly larger 
uncertainties and can be here ignored. 

We here compare these PAMELA preliminary data with the univocal predictions of the 
automatically stable Minimal Dark Matter candidate: the fermion quintuplet. Its contribution 
was predicted in the present paper in figures 5 and 6. In fig. 7 we sum it to the expected 
astrophysical background, assuming the standard NFW halo profile and the 'medium' model 
for positron diffusion in our galaxy. 

We see that the anomaly in positrons can be well reproduced; thanks to the Sommerfeld 
enhancement of the DM DM annihilation cross section in the non-relativistic limit, a modest 
boost factor is needed, B = 30 for the central value of M — 9.6 TeV suggested by the cosmo- 
logical DM abundance, and 3 < B < 100 in the 3a range of M compatible with cosmology, 
9.2 TeV < M < 10.2 TeV. 

In fig 7b we show the corresponding effect in the anti-proton flux. Although other possibili- 
ties exist, for simplicity we assumed the same energy-independent boost factor as for positrons. 
In such anomaly in anti-protons is expected at higher energies, compatibly with 

present PAMELA data. In the future, PAMELA should provide experimental data up to 
maybe 200 GeV for positrons, and 100 GeV for anti-protons. Higher energy will be reached by 
the future AMS experiment. These measurements will test the predictions illustrated in fig. 7. 
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To illustrate the discriminatory power of such future data, we point out an alternative 
interpretation of the PAMELA anomaly. In order to minimize the needed boost factor, we 
consider an hypothetical ad hoc DM candidate with mass M just above 60 GeV that annihilates 
directly into e~e + with 100% branching ratio. Assuming negligible co-annihilations and that 
DM DM annihilations are s-wave dominated, the cosmological abundance is reproduced for 
(av) « 3-10~ 26 cm 3 /sec. Inserting dN e + /dE e + = 5(E e +-M) in eq. (12), and taking into account 
positron energy loss in our galaxy as in eq. (13), we find that a boost factor B w 10 gives a 
positron fraction that somewhat resembles the PAMELA anomaly. If this other interpretation 
is true, PAMELA will see that the positron excess abruptly terminates at E e + > M. 
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